Striped morphologies induced by magnetic impurities in d-wave 

superconductors 

Xian-Jun ZucQ 

Department of Physics, Zhejiang University of Technology, 
Hangzhou 310023, Peoples Republic of China 
(Dated: November 22, 2010) 

Abstract 

We study striped morphologies induced by magnetic impurities in d-wave superconductors 
(DSCs) near optimal doping by self-consistently solving the Bogoliubov-de Gennes equations based 
on the t — t' — U — V model. For the single impurity case, it is found that the stable ground state 
is a modulated checkerboard pattern. For the two-impurity case, the stripe-like structures in order 
parameters are induced due to the impurity-pinning effect. The modulations of DSC and charge 
orders share the same period of four lattice constants (4a), which is half the period of modula- 
tions in the coexisting spin order. Interestingly, when three or more impurities are inserted, the 
impurities could induce more complex striped morphologies due to quantum interference. Further 
experiments of magnetic impurity substitution in DSCs are expected to check these results. 
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I. INTRODUCTION 



Recent studies of spatial inhomogeneous phases with stripe or checkerboard modula- 
tions in copper oxide-based compounds has received much attention. Neutron scatter- 
ing (NS)experiments reveal the existence of incommensurate magnetic peaks in cuprates, 
which has led to discussions of a stripe phase Q-Q. Meanwhile, scanning tunneling mi- 
croscop y (S TM) experiments observed checkerboard-like charge-density wave (CDW) modu- 



lations 



12 



14j . which break the square symmetry of the underlying lattice. It was proposed 



that the effect can be understood in terms of quasiparticle interference (QPI) due to scat 



tering on impurities and other inhomogeneities 15Ml8| within the so-called 'octet' model, 



or in terms of static or fluctuating stripes [191 ] . Kohsaka et al. suggest that both dispersive 
and non-dispersive modulated patterns originate from different regions in momentum and 
energy space [ijj]. The well-defined states in the nodal region are responsible for the low 
energy QPI structure, whereas the ill-defined quasiparticle states in the antinodal regions 
are responsible for the non-dispersive charge order above some energy scale. On the other 
hand, zero-field experiments on the well ordered YBCO samples show no evidence of static 
broken translational symmetry in nuclear magnetic resonance or NS experiments, which 
leads to the proposal that the checkerboard patterns observed by STM is caused by the 



local disorders 



20] . In brief, the physics that determines the ultimate patterns of competing 



orders is subtle. For instance, there are striped structures interpreted in terms of actual or 

n 

incipient orders [19J, or attributed to be induced by lattice distortion or impurity-pinning 



effect 



20 



21] . To date, the issue of the nature of spatial inhomogeneous phases in cuprates 



is still controversial. 

Impurity effects have been proven to be a valuable probe to explore the fundamental 
properties of cuprates since it often provides important insights into the underlying physics. 
Experimentally, STM detection of the modulation of the local density of states by impurities 
was used to great advantage to probe the nature of quasi-particle states of DSCs both in 
the superconducting and pseudo-gap phases [22 24 j . Another interesting issue is interfer- 
ence patterns and pinning stripes by impurities in DSCs. In general, impurity substitution 
causes additional slowing down of spin fluctuations and pinning of the stripes, leading to 
the formation of a static charge or spin order. Zhu et al. have proposed an explanation 
to the checkerboard pattern around a impurity or vortex cores based on an effective mean- 
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field t — U — V model [15j. For the case without an applied magnetic field, the disorder 
can produce a similar pinning effect of the fluctuating stripes. Andersen et al. investigate 
disorder-induced freezing of incommensurate spin fluctuations, which agrees qualitatively 
with experimental observations |25] . In the present work, we study the striped morphologies 
induced by magnetic impurities in d-wave superconductors (DSCs) near optimal doping. 
Including the competitions and coexistence among the DSC, spin-density wave (SDW), and 
CDW orders, the system is explored by self-consistently solving the Bogoliubov - de Gennes 
(BdG) equations. It is found that the local disorders can induce interesting phenomena, 
including checkerboard, stripe modulations, and other complex striped morphologies. We 
expect that these phenomena could be observed in the STM and NS experiments. 

II. THEORY AND CALCULATION DETAILS 

We start from the two-dimensional t — t' — U — V model, which consists of two parts, 
H = H Q + H imp . The Hamiltonian H and H imp describe the superconductor and magnetic 
impurities, respectively, which can be written as 

H Q = -5^j(cLc i<T + H.c.) 

ija 

ij ia 

H imp = J2h eff (i)(4 t c n - cj ; c 4 ). (1) 

i 

Here Cj CT annihilates an electron of spin a at the ith site. The hopping integral takes t 
between nearest neighbor (NN) pairs, and t' between next-nearest neighbor (NNN) pairs. U 
is the on-site Coulomb repulsion interaction, fi is the chemical potential, which is determined 
self-consistently in the calculation. The local effective field h e ff(i) is introduced to model 
the exchange coupling between conducting electrons and the impurity spin, where we have 
treated the Kondo impurity spin as a Ising-like one. Some similar model had been employed 
to study the effects of magnetic impurities on cuprate superconductors [26l . |27| ] , which can 



qualitatively explain the observed impurity states well 23]. Therefore, we employed the 
above model in this work. The self-consistent mean-field parameters are given by = 
Y^ a < c\ a c ia >, the magnetization = (l/2)(< cLc^ > — < c^c^ >), and the DSC order 
parameter Ay = (V/2) < c^Cj^ — c^c^ > with V the phenomenological pairing interaction. 



The Hamiltonian H can be diagonalized by solving the following BdG equations, 



^ -//;,. 



*j = E%, (2) 



where the quasiparticle wave function = (u^, Vi^) T . The spin-dependent single-particle 
Hamiltonian reads Hij a = — t8i +T j— i'^i+r'j + Ej cr ^e//(«)^,i m + ^i,a-— Here the sub- 
scripts r and r' denote the unit vector directing along four NN and NNN bonds respectively, 
and i m is the position of the impurity site. The self-consistent parameters are given by = 
£J< t | 2 /(£„), n iX = £ n KI 2 [l - f(En)), and A y = \ £>? t ^ + v^]tanh(^), 
where f(E) = 1/(1 + e l3E ) is the Fermi-Dirac distribution function. Hereafter, the length is 
measured in units of the lattice constant a, and the energy in units of t. We set U = 2.5 
and t' = —0.2 in this paper. The pairing interaction is chosen as V = 1.0 to guarantee 
that the superconducting order A — 0.08t, comparable with the observed T c in cuprate 
superconductors. We study the impurity effects on the electronic states of DSCs near op- 
timal doping with the filling factor rij = ^ io . c\ a c i(T / (N x N y ) = 0.83 (i.e., the hole doping 
x = 0.17), where N x , N y are the linear dimension of the unit cell. The BdG equations 
are solved self-consistently for a square lattice of 32 x 32 sites, and the periodic boundary 
conditions are adopted. The numerical calculation is performed at a very low temperature, 
T = 10 _5 K, to extract the low-energy physics. The local effective field is taken to be h e ff at 
the impurity site and zero otherwise. The DSC order parameter at the ith site is defined as 
Aj = (A lii+eT + A^i-c - A i)i+eH - A M _ e J/4, and the spin order parameter is Mi = (-1)^. 



III. NUMERICAL RESULTS AND DISCUSSION 

In FIG. [U we plot the spatial distributions of the DSC, spin, and charge orders around 
the impurity site. One can see that all the three orders display checkerboard modulations 
around the magnetic impurity. Similar to the nonmagnetic impurity case, a SDW checker- 



board pattern is observed around the magnetic impurity 



15 



which coincides with the 



n 

NS data However, it is noteworthy that the checkerboard pattern of the DSC order can 
also be induced by the magnetic impurity, and a weak associated CDW pattern is observed, 
which is different from the nonmagnetic impurity case. Moreover, the modulated DSC and 
CDW orders share the same periodicity. Overall, the DSC order is strongly suppressed at 
the impurity site while the amplitudes of the CDW and SDW orders reach global maxima. 



This is the common feature of orders around the magnetic impurity despite various param- 
eters. Therefore, in view of these phenomena, we clearly see the relationship of competition 
and coexistence between antiferromagnetic and DSC orderings. In comparison to the non- 
one finds that in both cases the DSC orders are suppressed 



magnetic impurity case 
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around the impurity site. However, for the CDW order, the opposite tendency is observed 
due to the fact that the nonmagnetic impurity is repulsive. 

Next we consider the many-impurity case to show the striped morphologies which could 
be induced by magnetic impurities in DSCs. As shown in FIG. [2j upon inserting magnetic 
impurities in the same direction, a local stripe-like structure can be induced, which is analo- 



gous to the nonmagnetic impurity case [281 ] . A second important feature is that the average 
separation between neighboring stripes is not expected to change upon impurity-doping. 
The impurity-pinned SDW stripe also show a modulation with periodicity 8a, which coex- 
ists with the DSC order. In addition, similar to the single impurity case, the modulated 
DSC and CDW stripes also share the same periodicity. For the three-impurity case, one 
notes that actually there are local ferromagnetic moments formed around impurities. If the 
three impurities are placed so as to form a right-angled shape, one finds that the striped 
structures also show a right-angled shape due to the impurity-pinning and quantum inter- 
ference effects (see Fig. 3). The impurity-pinned stripes along x or y direction are balanced 
by quantum interference effects so that structures with a right-angled shape are favorable 
on energy. Meanwhile, the alternately dark and bright features due to quantum interference 
effect are still obvious. Away from the impurities, the stripes eventually evolve to checker- 
board structures due to quasiparticle interference, and tend to disappear on farther sites. 
This feature can be seen more clearly from the DSC order parameter [see Fig. 3 (a)]. 

Moreover, we also calculate the cases with four or five impurities. As shown in Fig. HI 
when inserting four or five magnetic impurities symmetrically around the center site, there 
are symmetrical impurity-pinned stripe structures formed around impurities. Interestingly, 
these stripe structures are somewhat analogous to "quantum corrals", structures built from 
adatoms using atomic manipulation processes, which have been observed on a metal sur- 



face 



29 



301 ] . Here, these quantum-corral-like structures in order parameters are produced 



due to impurity-pinning and quantum interference effects in DSCs. This feature is espe- 
cially clear in the spin order, in which the alternately dark and bright rectangular corrals 
surround the impurities near the center. Similar to the three-impurity case, the stripes even- 
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tually evolve to checkerboard structures in the sites away from the impurities and disappear 
on farther sites. 

IV. CONCLUSION 

To conclude, we have investigated the modulations around magnetic impurities of the 
DSC, spin, and charge orders in DSCs. It is found that magnetic impurities could induce 
rich and interesting striped morphologies. For the single impurity case, the checkerboard 
modulations are observed in all three orders. The weak associated CDW pattern is dif- 
ferent from the nonmagnetic impurity case. When more magnetic impurities are inserted, 
more complex modulated structures could be induced, including rectilinear and right-angled 
stripes and quantum-corral-like structures due to impurity-pinning and quantum interference 
effects. The latter structures were only reported on metal surfaces previously. We expect 
that these phenomena could be observed in the high resolution STM and NS measurements 
on DSCs. 
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FIG. 1: (Color online) The surface plots of orders around the magnetic impurity on a unit cell of 
size 32 x 32 sites, (a), (c), and (e) are the spatial distributions of the DSC, spin and charge orders, 
(b), (d), and (f) are their contour plots. 




FIG. 2: (Color online) The spatial contour plots of orders by inserting two or three magnetic 
impurities in the same direction, (a), (c), and (e) are the contour plots of the DSC, spin and 
charge orders with two impurities, while (b), (d), and (f) are the three-impurity case. 
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FIG. 3: (Color online) The same plot as Fig. 1 but with three impurities forming a right-angled 
shape. 
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FIG. 4: (Color online) The spatial contour plots of orders by inserting four or five magnetic 
impurities symmetrically around the center site, (a), (c), and (e) are the contour plots of the DSC, 
spin and charge orders with four impurities, while (b), (d), and (f) are the five-impurity case. 
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